Competing mechanisms for step meandering in unstable growth 
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The meander instability of a vicinal surface growing under step flow conditions is studied within 
a solid-on-solid model. In the absence of edge diffusion the selected meander wavelength agrees 
quantitatively with the continuum linear stability analysis of Bales and Zangwill [Phys. Rev. B 
41, 4400 (1990)]. In the presence of edge diffusion a local instability mechanism related to kink 
rounding barriers dominates, and the meander wavelength is set by one-dimensional nucleation. 
The long-time behavior of the meander amplitude differs in the two cases, and disagrees with the 
predictions of a nonlinear step evolution equation [O. Pierre-Louis et al., Phys. Rev. Lett. 80, 4221 
(1998)]. The variation of the meander wavelength with the deposition flux and with the activation 
barriers for step adatom detachment and step crossing (the Ehrlich-Schwoebel barrier) is studied in 
detail. The interpretation of recent experiments on surfaces vicinal to Cu(100) [T. Maroutian et al., 
Phys. Rev. B 64, 165401 (2001)] in the light of our results yields an estimate for the kink rounding 
barrier at the close packed steps. 



I. INTRODUCTION 

It has been shown in several experiments and com- 
puter simulations that during epitaxial growth on a vici- 
nal crystal surface, straight steps amomstable against the 
formation of an in-phase meanderEIEI. The phenomenon 
was first predicted theoretically by Bales and Zangwill 
(BZ) within a continuum theorytl According to BZ, step 
meandering is caused by, an energy barrier, the Ehrlich- 
Schwoebel (ES) barrieiu, which suppresses the attach- 
ment of surface atoms to the step from the terrace above. 
The preferential attachment from below implies that pro- 
trusions in the step are amplified, leading to a linear dif- 
fusional instability. 

Recent experimental measurements xif. the meander 
wavelength on vicinal copper surfacestm show signifi- 
cant disagreement with the predictions of the BZ the- 
ory. This -has. lead to a search for alternative sources of 
instabilityLjfjiJ. The most prominent alternative mecha- 
nism is the kink Ehrlich-Schwoebel effect (KESE), which 
invokes a kinetic barrier that prevents atoms-diffusing 
along step edges to cross corners or kinksEiJOo. In 
close ariaJpflw-|to the ES instability of singular crystal 
surfacegl3li3'L5l, this induces a destabilizing mass current 
along the step. 

In this article we present kinetic Monte Carlo (KMC) 
simulations of step meandering which display both types 
of instability within a single model. The relative impor- 
tance of the KESE vs. the BZ instability can be tuned 
through the facility of step edge diffusion. By explicitly 
relating the parameters of the KMC model to those of 
the continuum theory, we show that the meander wave- 
length can be quantitatively predicted in both instabil- 
ity regimes. Our study thus proves the feasibility of ex- 
tracting kinetic barriers from experimental meander mor- 
phologies. The simulations also provide some insight into 



the long-time behavior of the meander amplitude which 
can be compaij cd to t he predictions of nonlinear contin- 
uum equationsMliZllla. 

The model employed in our work is described in the 
next section. Section III summarizes the predictions of 



continuum theory for the meander wavelength, and ex- 
plains how the material parameters of the continuum de- 
scription are determined for the KMC model. The sim- 
ulation results are presented in Section |y|. We provide 
some discussion of the applicability of the KESE scenario 
to the experiments^ on Cu(100) in Section |v|, and con- 



clude in Section VI 



II. MODEL 

We employ a standard solid-on-solid modclE^I, in which 
the crystal surface is represented by columns of integer 
height h r on a square lattice of substrate sites r. The 
elementary processes are the deposition of atoms at rate 
F and the hopping of adatoms to nearest neighbor sites 
with a rate 



r exp(-£: a /£; B T). 



(1) 



The activation barrier E a depends on the local configu- 
ration through 

E a = Es + n l E n + (rii - n/)9(rt l - n f )E B B (2) 
+ (m, - m/)9(r7ij - m f )E E s, 

where E§ is the energy barrier for diffusion on a flat ter- 
race, E n is the contribution of a nearest neighbor bond 
to the energy barrier, Ebb is an additional energy cost 
for bond breaking and £"es is the ES barrier; nt denotes 
the number of in-plane nearest neighbors before the hop 
and rif after the hop, while mi,mf are the number of 
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next-nearest neighbors in the planes beneath and above 
before (rrii) and after (rrif) the hop. The Heaviside func- 
tion 9(a;) = 1 if x > and otherwise. 

The implementation of the ES barrier through the 
change in the number of out-of-plane next nearest jpgigh- 
bors has been used in several earlier growth studieso'EllEiJ, 
and the additional bond breaking energy. _Ebb was intro- 
duced in the context of ion sputteringEj. The rates de- 
fined by (Q) and (Q) satisfy detailed balance with respect 
to the Hamiltonian 
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The sum runs over all nearest neighbor pairs, and 
Ek = ^ (En + E B b ) 



(4) 



is the energy per unit length of a single height step run- 
ning along one of the lattice axes; for the SOS model, this 
is also the kink energy. The detailed balance condition 
is easily checked by noting that the rates can be writ- 
ten as a product of Arrhenius (term proportional to E n 
in @) and Metropolis (terms proportional to -Ebb and 
Ees) dynamicsEil, each of which fulfil detailed balance 
with respect to part of the Hamiltonian (||). 

Setting Ebb = the model, called hereafter model I, 
does not include diffusion along the step edges, because 
the hopping rate along the step is equal to the rate of 
detachment from the step. Edge diffusion is facilitated 
compared to detachment if E-qb > (called hereafter 
model II). Model II also contains a kink ES barrier, since 
atoms cannot round corners without detaching from the 
edge. 

The simulations were carried out on rectangular lat- 
tices with periodic boundary conditions in the step di- 
rection and helical boundary conditions in the direction 
of the vicinality. The initial step spacing was typically 
I = 6 (exceptions are noted in the figure captions). For 
both models the values for the activation barriers were 
set to E s = 0.35 eV, E n = 0.21 eV and E ES = 0.15 
eV, Ebb — for model I and for model II we put 
Ebb = E n = 0.21 cV. The temperature was T = 375 K, 
the diffusion prefactor r = 2xl0 11 s~ 1 (r = 4xl0 12 s~ 1 ) 
for model I (model II) and the deposition flux was varied 
in the range F = 0.005 — 1.0 ML/s. These choices were 
motivated mainly by our desire to study the variation of 
the meander wavelength over a range of control parame- 
ters without being strongly affected by finite size effects 
and limited computer time. In particular, the different 
diffusion prefactors for the two models were chosen only 
to be able to study both models in the same range of flux 
and temperature, and do not carry any physical signifi- 
cance. 

Typical configurations generated in the simulations are 
shown in FigsJ^J| and || For both models the initially 
straight steps form an in-phase meander pattern with a 



characteristic wavelength. The selected wavelength re- 
mains constant during growth. The dependence of the 
meander wavelength on the model parameters is the main 
focus of the following discussion. Some aspects of the 
tem poral evolution of the pattern will be addressed in 

SectJivrj. 



III. RELEVANT LENGTH SCALES 

Before turning to the quantitative analysis of the sim- 
ulations, we summarize the available theoretical predic- 
tions for the length scale of the meander instability. 

A. The Bales-Zangwill instability 

The BZ analysis proceeds by solving a diffusion equa- 
tion for the adatom concentration with boundary con- 
ditions given by the attachment-detachment—kinetics at 
the steps. It predicts_.an in-phase meanderEj with the 
dominant wavelengthlia 
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Here ft denotes the atomic area and a is the lattice spac- 
ing. All kinetic and thermodynamic parameters entering 
(||) can be expressed in terms of the microscopic rates 
and energies of the KMC model. The diffusion coefficient 
on the terrace reads D = c e .exp(— Es/kgT), and the 
step edge stiffness is given byE4 T = 2asmh 2 (E K /2k B T). 
To calculate the equilibrium adatom concentration cjL 
and the mobility a s t along the step edge, it is useful 
to consider the step as a one-dimensional (ID) SOS in- 
terface in equilibrium with the adatom gas on the ter- 
race. For model I the transition rates are of Arrhe- 
nius form, i.e. dependent only on the configuration be- 
fore the iump, and exact results for c° eq and a st can 
be found!! We obtain c° eq = exp(-2E n /k B T) and 
a st = (D/2)cxp(-2E n /k B T). 

For model II the transition rates-are no longer of Ar- 
rhenius type, and the exact resultsa cannot be directly 
applied. Within linear fluctuation theory, the step edge 
mobility has been estimated as a st = a 2 /t£, where 3U, 
is the characteristic time for detachment from a kinko. 
For our model II this yields 



a st « £exp [-{2E n + E B B)/k B T} 



(6) 



Strictly speaking, (J^) has to be modified in the presence 
of kink ES barriers which reduce mass transport along 
the edge. A detailed analysisEH shows that an additional 
KESE barrier Ekes becomes relevant when Ekes > Ek- 
In our simulations Ekes ~ Ebb < Ek always, so (|^) 
suffices. The terrace adatom concentration is generally 
given by the expression c° eq = exp(— AB/fcgT), where 
AE = AEk is the formatioftjenergy for moving an adatom 
from a kink to the terraceEa. 



FIG. 1: Evolution of a vicinal surface with initially straight 
steps under the dynamics of model I. Snapshots show a 120 x 
500 piece of a larger system ( 120 x 1000 with 20 steps) after 
deposition of 50, 100 and 200 monolayers. The deposition 
flux was F = 0.2 monolayers per second (ML/s), and other 
parameters as described in the text. 



FIG. 2: Evolution under the dynamics of model II. Snapshots 
show a 120 x 500 piece of a larger system (120 x 800 with 
20 steps) after deposition of 5, 20 and 60 monolayers. The 
deposition flux was F = 0.01 ML/s and other parameters as 
described in the text. 



The strength of the ES effect is contained in the pa- 
rameter.^ = (Z_ — + I- + l+), where the length 
scalesEflEj l± = D/k± are inversely proportional to the 
attachment rates for adatoms approaching the step from 
below (fc + ) or above (fc_) (Z_ is also knownE_3 as the ES 
length Ies)- For the KMC model we estimate 



k+=D , w Dc X p(-EEs/k B T). 



(7) 



In fact the rate of attachment from the upper terrace 
depends on the microscopic configuration of the step. 
A direct test of the relation (uv is presented below in 



Sect. IV C 



The expression used above for the edge mobility a s t is 
valid only on length scales longer than the kink distance 
Ik ~ (1/2) exp(Ex/kBT). On shorter scales edge diffu- 
sion is much more efficient, and the expression for the 



characteristic meander wavelength must be replaced byEI 

Xbz = 2 1/4 yfh^ xTz (8) 
when edge diffusion dominates and Ik S> ^bz- 

B. The KESE instability 

Step meandering due to the KESE can be discussed0 
in ana] nay, to the ES instability of a one-dimensional 
surfacell3ll3. The characteristic wavelength depends on 
the one-dimensional nucleation length lo, which is de- 
fined as the average distance between two dimers that 
are nucleated on a straight step at the beginning of de- 
position, and the kink ES length Zkes ~ expfEKEs/^sTl- 
For a strong KESE, in the sense that Id Zkes> the ini- 
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tial meander wavelength is Id, r^We for a weak KESE 
(In > Jkes) it is of the order ofHLa 

A„ « (£ D /t K Es) 1/2 iD- (9) 
From one-dimensional nucleation theory the expression 
l D = {UDst/Fl) 1 ^ (10) 

for the nucleation length can be derivedS; here Fl is the 
flux of terrace atoms onto the step. In contrast to the 
mobility a s t, the edge diffusion coefficient D st refers to 
the motion of an edge atom along a straight step without 
kinks, and is given by D st — D exp(— -E^/feaT) in the 
SOS model. 

It should be noted that (^,|l^) apply to a ID surface in 
the absence of desorption; for a step this translates into 
neglecting the detachment from the step. This approxi- 
mation clearly breaks down as the bond breaking barrier 
Ebb approaches zero. Including the detachment from 
the step introduces a new length scale into the problem, 
namely the diffusion length x s — \J D s tT. Here r is the 
time an adatom diffuses along a straight step before being 
"evaporated" to the terrace, and x s the distance the atom 
travels along the step in time interval r. A lower bound 
for r is given by the detachment rate from the straight 
step r -1 = D st cxp(— E-qb /ksT). In reality detached 
atoms have a high probability of re-attachment and thus 
the real evaporation time r is longer. A lower bound for 
the diffusion length then reads x s — exp(i?BB /2ksT). 
Repeating the calculation Qf-,the nucleation length for a 
ID surface with desorptionEH, one arrives at an expres- 
sion for the nucleation length in the limit x s <C Id, which 
reads 

1/2 

Id=[^) -, (11) 



The same result can be obtained .also by applying the 
scaling arguments of Jensen et alEB in one dimension. 
For general values of x s , Id is given as the solution of a 
transcendental equation. A simple expression interpolat- 
ing between the two cases ( |To| ) and (|ll]) reads 



In = 



Dst 
Fl 



1/4 
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(12) 



In the case of a weak KESE, repeating the calculation of 
Politi and VillairJij with desorption, one findsED that for 
x s <c Ikes ^ Id the most unstable wavelength is of the 
order 



A,, 



l D / x s 



(13) 



which replaces 



IV. SIMULATION RESULTS 

A. Meander mechanisms 

Figure |] shows the meander wavelength as a function 
of deposition flux obtained from our KMC simulations. 
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FIG. 3: Step configurations for model I (upper panel) and 
model II (lower panel) after deposition of 75 ML at a deposi- 
tion flux F = 0.5 ML/s (model I) and F = 0.01 ML/s (model 
II). The figures show part of a 800 x 120 lattice with 20 steps. 



The wavelength was extracted from the profiles directly 
by measuring the distance between subsequent minima 
in a single profile. The error bars refer to the variation 
of the wavelength within a single profile. For model I 
the wavelength is found to scale as ^-0.47:1:0.06 ^ j n q ua i_ 
itative agreement with both the BZ-length Xbz and the 
nucleation length Id in the detachment-dominated limit 
(Eq.(pl]); for model I x s w 1 <C Id)- Quantitatively the 
results are found to agree with Xbz for the parameters 
used in the simulations, while the nucleation length Id 
is smaller by approximately a factor 1/2. A more con- 
vincing piece of evidence for the Bales-Zangwill mecha- 
nism is the dependence of the meander wavelength on the 
Ehrlich-Sc hwoe bel barrier E-^s, which is discussed below 
in Section IV C . The nucleation length is obviously inde- 
pendent of Efts- 

For model II the meander wavelength scales as 
F~°- 28±0 - 05 , which disagrees with the BZ theory but is 
consistent both with the modified BZ length (^), and 
with the nucleation length ( |Io| ) in the absence of detach- 
ment. However, Eq.(^|) predicts a prefactor that is one 
order of magnitude too large. This is not surprising, since 
(||) was derived for steps close to thermal equilibriumLl; 
under growth conditions the kink density is much larger 
than its equilibrium value. 

Including detachment and using the lower bound x s — 
exp(£ , Bs/2fcsT) as an approximation for the diffusion 
length yields an upper bound for the nucleation length. 
Expressions ( ]l0| ) and ( |l2"| ) thus give lower and upper 
bounds for the nucleation length, which differ approxi- 
mately by a factor of 2. The lower bound ( |l0| ) is seen 
to quantitatively describe the simulation data for model 
II, which shows that for model II the detachment may 
be neglected. This confirms that one-dimensional nu- 
cleation is the relevant wavelength selection mechanism 
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FIG. 4: Meander wavelength as a function of flux for model 
I (circles) and model II (diamonds). Each symbol represents 
a single run on a lattice of size 1000 x 30 with 5 steps. The 
error bars refer to the variation of the wavelength within the 
configuration. For some fluxes results for a lattice of size 30 x 
1200 have been included also. The simulations were run until 
the meander wavelength was clearly visible. The BZ-length 
(^) is plotted as a full line for model I and a short-dashed line 
for model II. The long-dashed line is the nucleation length 
@ for model II. 



under conditions of facile edge diffusion, in accordance 
with the conclusions from previous experimentalEl and 
simulationaloB work on surfaces vicinal to Cu(100). 

The meander wavelength by itself does however not 
uniquely specify the instability mechanism. The BZ the- 
ory predicts a band of unstable wavelengths extending 
over the interval (X m i n ,oo), where \ m in = Xbz/V% and 
Xbz, as given by (|j), is the wavelength of perturba- 
tions with the maximal growth rate. Numerical studies 
of a nonlinear evolution equation for the in-phase me- 
ander show that initial wavelengths between X m in and 
Xmax ~ 3Xbz are preserved during growthliH. Thus de- 
viations from the BZ prediction (^J) can be attributed 
partly to a wavelength different from Xbz which domi- 
nates the spectrum of initial perturbations. In this con- 
text it is important to note that, for small fluxes, the 
wavelength measured for model II is smaller than the 
minimal unstable wavelength X m i n of BZ theory; this can 
be seen by shifting the short-dashed line in Fig.[|, which 
represents Xbz for model II, downward by a factor of \/2. 
This proves that for model II an instability mechanism 
different from the BZ mechanism - the KESE - is the 
cause of the meander. 



B. Crossover between the two mechanisms 

Which of the two meander mechanisms is operative 
depends on the importance of step edge diffusion and on 
the strength of the KESE barrier. In our SOS model 
both are controlled by the bond breaking barrier Ebb ■ 
Thus decreasing Ebb should lead to a crossover from the 
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FIG. 5: The observed wavelengths (circles) for various values 
of the bond breaking barrier Ebb ■ The full line is the nucle- 
ation length (|l2|), the dashed line the BZ-length (^) and the 
dot-dashed line the meander wavelength (|l^) for weak KESE. 
The crossover between the two meander mechanisms occurs 
around Ebb ~ 0.05 eV. Each point is an average over 5 in- 
dependent runs on a 250. .1500 x 30 lattice with 3 steps (step 
spacing I = 10). 



KESE instability to the BZ instability. 

Simulation results for various values of Ebb are shown 
in Fig. |5|. The wavelength was determined by count- 
ing the number of minima in the step profile on a finite 
sample. Reported wavelengths are averages over 5 inde- 
pendent runs and the error bars are the standard devi- 
ations. Since the bond breaking barrier enters the kink 
energy through (||), the predicted length scales (||) and 
( |l2| ) both increase with decreasing -Ebb m a qualitatively 
similar manner. In addition, the decrease of the KESE 
barrier implies a transition from the strong KESE to the 
weak KESE regime. The numerical data clearly show 
that for Ebb > 0.10 eV the meander wavelength is set by 
the nucleation length Irj, while for Ebb < 0.05 eV the 
simulations are consistent only with the BZ length, be- 
cause the nucleation length is considerably smaller, and 
the weak KESE length ( |l3|) is much larger than the ac- 
tual meander wavelength. 



C. Variation of the ES barrier 

As a further test of the BZ prediction (|5|) we have mea- 
sured the meander wavelength in model I for different 
values of the ES barrier. At T = 375 K and with the 
barrier Ees = 0.15 eV used above, the Schwoebel length 
Z_ w 104 > I and l + = 1, so that fg « 1 in Eq. (||). 
To study the effect of a finite ES. length we carried out 
simulations with a single terraceEd of width I = 15, vary- 
ing the ES barrier between E E s = 0.06 eV and 0.20 eV 
(/_ = 6 — 487), while keeping the other barriers at the 
values given previously. Again, the average wavelengths 
were calculated from 5 independent runs. Figure ^ shows 
that the dependence of the meander wavelength on the 
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FIG. 6: Meander wavelength as a function of the ES length 
Ies = D/k- for model I at F = 0.25 ML/s. Each symbol is 
an average over 5 runs on a lattice of size 1250 x 15 containing 
a single step (step spacing I = 15). The dashed line shows the 
BZ prediction (Bp. 
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FIG. 7: Average step width w as a function of time for model 

I (upper panel, F = 0.1(O), 0.4(D), 0.8(o) ML/s) and model 

II (lower panel, F = 0.05(*), 0.2(o), 0.4(D), 0.8(A) ML/s), 
from single runs on a lattice of size 500 x 30 with 5 steps. 



ES barrier is qualitatively described by BZ theory, but 
Eq. (||) is not quantitatively accurate. The true ES length 
appears to be smaller than that given in (^) by about a 
factor 1/4, which corresponds to a reduction of the ES 
barrier by 0.04 eV. This cannot be a simple effect of step 
roughness, since the implementation of the ES barrier 
used in the present work in fact implies that the ap- 
proach to a kinked step involves a higher barrier than 
to a straight step. This issue deserves further considera- 
tion. 



D. Temporal evolution 

The two models can be distinguished with regard to 
the dynamics of meander formation. For model I the 
steps meander in-phase from the beginning, whereas for 
model II the meander starts with random phase shifts be- 
tween the steps (compare the early time configurations in 
Figures |l| and ) . Later on, as the meander amplitude 
grows and the step-step interaction through the diffu- 
sion field becomes effective, the correlations between the 
steps grow, leading asymptotically to an in-phase step 
train. In the light of the two different mechanisms these 
observations are easy to understand: Since for model I 
the meander is due to the BZ instability, the step train 
is expected to be in-phase iijom the beginning, this being 
the fastest growing modeEj; for model II the meander 
starts independently at each step since the meander is 
due to the local adatom dynamics at the step edge. 

We turn next to the time evolution of the meander 
amplitude, which also differs for the two models (Figj^). 
We consider the step width w defined by 

L 

w 2 = L-^[y(x)} 2 , (14) 

x=l 



where the cc-coordinate is directed along the step, L is 
the step length and y(x) is the step position relative to 
its mean. For model II, w(t) increases verv-jrapidly at 
the beginning, but asymptotically saturatesE^I. In con- 
trast, for model I the step width was found to increase 
linearly in time for the longest times we could accessEa. 
Both observations are at variance with the predictions 
of a nonlmeax-step evolution equation for the in-phase 
meanderEjilZl'tj, which predicts that w ~ \ft. 

Our results also contradict earlier Monte-Carlo 
simulationllB, in which the steps were described by sin- 
gle valued functions, thus prohibiting step overhangs and 
voids, and nucleation on the terraces was not allowed 
for. In these simulations a regime with w ~ y/t was ob- 
served at long times. A direct comparison of the two 
sets of simulations is not possible, however, because of 
the rather different choice of parameters. In the earlier 
workllS, the concentration of adatoms was set to a value 
c Q eq = 0.119, which is several orders of magnitude larger 
than in our simulations. In fact, at such high adatom 
concentrations step flow is hardly possible, if the nucle- 
ation on the terraces is not artificially suppressed. For 
this reason we cannotr^eproduce the conditions used in 
the earlier simulationsEj. 

The linear amplitude growth in model I suggests that 
minima and maxima of the step meander move at differ- 
ent velocities. This may be due to diffusional screening 
which prevents adatoms to reach the narrow fjords sep- 
arating the cells in the upper panel of Fig|^. In model 
II the fjords are wider, because adatoms are able to fill 
them up through edge diffusion. The saturation of the 
amplitude in model II may be related to the stabilizing 
character of the KESE step current at large slopeai3'E3. 
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V. RELATION TO EXPERIMENTS ON CU(IOO) 

In this section we briefly comment on the relevance 
of our wqik for the experimentally observed meander 
instability^ on surfaces vicinal to Cu(100). In these ex- 
periments two different vicinal surfaces were considered, 
which consist of dense packed (110) steps, and open (100) 
steps, respectively. 

For the dense packed (110) steps, edge diffusion is 
much more facile than detachment, so the scenario of 
our model II should apply. Indeed, using the expression 
© for the one-dimensional nucleation length to inter- 
pret the experimentally measured activation energy E a = 
0.092 eV of the meander wavelength, A ~ e ~ E a/k B T ^ one 
obtains an energy barrier E st = 4E a — 0.37 eV for diffu- 
sion along a straight edge, which is consistent with the 
estimate E st — 0.45 ± 0.08 eV derived, from the analysis 
of time-dependent STM observations^. Also the flux de- 
pendence of (|l(J) as E -1 / 4 agrees with the experimental 
power law exponent of —0.21 ± 0.08. 

What is missing to complete the picture is some di- 
rect experimental evidence for a (strong) KESE at the 
(110) stcpEi Here we want to point out that indirect 
evidence for a kink ES barrier follows from a comparison 
of the growth experimentscl with stpc fluctuation mea- 
surements. Using the accepted valueo £k ~ 0.13 eV for 
the kink energy, the measurement f& the prefactor of the 
temporal step correlation functional yields the estimate 
E a 0.91 eV for the activation energy of the step edge 
mobility o~ s t- It was mentioned above in Section III that, 



in the absence of a strong KESE, this can be identified 
with the energy barrier -Edet for the detachment of a step 
adatom from a kink. In a simple bond counting approxij 
mation (which is supported bjL£ffective medium (EMT)Ea 
and embedded atom (EAM)Ej calculations) the detach- 
ment barrier is given by -Edct ~ E st + 2£k- Using the 
value of E st determined from the meander wavelength of 
the (110) step this yields E<\ t ~ 0.63 eV, which is much 
smaller than the step fluctuation estimate of E a . 

The discrepancy strongly suggests that the migration 
of atoms along the kinked step is suppressed by an addi. 
tional kink ES barrier Ekes- The quantitative analysisEZl 
shows that for Ekes > Ek, the activation energy for o~ s t 
is given by E a = Edct+ExES - Ek, which, using the num- 
bers quoted above, yields the estimate -Ekes tMM 1 
In agreement with semiempirical calculationsE3o, the 
additional kink barrier is found to be comparable to the 
barrier E st for diffusion along a straight step. 

The situation is rather different for the open (100) 
step, at which a meander instability with similar char- 
acteristics has been observed!! Maroutian et al. iden- 
tify the measured wavelength for the open step with the 
one-dimensional nucleation length (jl0|). This interpre- 
tation seems questionable for the following two reasons. 
First, since the open step can be viewed as being com- 
posed entirely of (110) kinks, the edge diffusion barrier 
E st is expected to be much larger than along the close 
packed stepc9, while the additional barrier for detach- 



ment from a kink (which-is related to the energy of (100) 
kinks) is much smallerEj. Under these conditions the 
initial growth of the step is not well described by the pic- 
ture of orLCadimcnsional nucleation. Second, symmetry 
argumental3ll3 show that the step edge current induced 
by the KESE should stabilize, rather than destabilize the 
open step. Thus the mechanism which could lead to an 
instability on the scale of the nucleation length is absent. 

This leaves the BZ instability as the most plausi- 
ble instability mechanismEil. However, neither the ac- 
tivation enecgyj— (as estimated from known energetic 
parameterst3c3'E2l) nor the E _1 / 2 -dependence predicted 
by (^) matches. J he experimentally determined mean- 
der wavelengtbtffl. In addition, in KMC simulations of 
Cu(100), the surface composed of (100) steps was found 
to be stableQ; possibly this reflects the stabilizing KESE. 
The experimentally observed instability of the open step 
thus poses a major puzzle at the moment. 



VI. CONCLUSION 

In summary, we have studied the unstable growth dy- 
namics of vicinal surfaces using an SOS model with two 
different choices of atomic processes. It was shown that 
in both cases the steps form an in-phase wave pattern. 
The formation mechanism and the wavelength of the pat- 
tern depends on whether the diffusion of adatoms along 
the step edges is significant or not. In the case where step 
edge diffusion is negligible the wavelength was found to 
be correctly predicted by the continuum theory of BZ, 
while in the case of facile edge diffusion (combined with 
an ES barrier at kinks) it is set by the one-dimensional 
nucleation length. 

In both cases good agreement between the KMC sim- 
ulations and the analytic predictions (^) and (|l0|) for the 
meander wavelength was achieved without any adjustable 
parameters. However, the asymptotic time evolution of 
the step profile disagrees with the prediction of an ef- 
fective step evolution equationliSEZllla. Whether this dis- 
agreement is due to different physics described by the 
continuum and discrete models or due to the approxi- 
mations made deriving the step equation remains to be 
clarified. 
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